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Abstract. We present the parameters of 891 stars, mostly clump giants, including atmospheric parameters, distances, absolute 
magnitudes, spatial velocities, galactic orbits and ages. One part of this sample consists of local giants, within 100 pc, with 
atmospheric parameters either estimated from our spectroscopic observations at high resolution and high signal-to-noise ratio, 
or retrieved from the literature. The other part of the sample includes 523 distant stars, spanning distances up to 1 kpc in the 
direction of the North Galactic Pole, for which we have estimated atmospheric parameters from high resolution but low signal- 
to-noise Echelle spectra. This new sample is kinematically unbiased, with well-defined boundaries in magnitude and colours. 
We revisit the basic properties of the Galactic thin disk as traced by clump giants. We find the metallicity distribution to be 
different from that of dwarfs, with less metal-rich stars. We find evidence for a vertical metallicity gradient of -0.31 dex kpc~' 
and for a transition at ~ 4-5 Gyr in both the metallicity and velocities. The age - metallicity relation (AMR), which exhibits a 
very low dispersion, increases smoothly from 10 to 4 Gyr, with a steeper increase for younger stars. The age-velocity relation 
(AVR) is characterized by the saturation of the V and W dispersions at 5 Gyr, and continuous heating in U. 

Key words. Stars: kinematics - Stars: fundamental parameters - Galaxy: disk - Galaxy: structure - Solar neighbourhood 



1. Introduction 

This paper is the continuation of previous papers (Soubiran 
et al. I20031 hereafter Paper I and Siebert et al. I2003I here- 
after Paper II) where we investigated the vertical distribution 
of disk stars with local and distant samples of clump giants. 
Our main result in Paper I was a new characterization of the 
thick disk, showing a rotational lag of -51 ± 5kms~' with 
respect to the Sun, a velocity ellipsoid of (cr^, <xy, ov) = 
(63 + 6,39 ± 4,39 ± 4)kms _1 , a mean metallicity of [Fe/H] 
=-0.48 + 0.05 and a high local normalization of 15 + 7%. We 
have also determined in Paper II the gravitational force perpen- 
dicular to the galactic plane and the mass density in the galactic 
plane (2 = 67M pc~ 2 ) and thickness of the disk (390^ P c )- 
We found no vertex deviation for old stars, consistent with an 
axisymetric Galaxy. After these two papers, we have enlarged 
and improved our samples in order to go further into the study 
of the local thin disk. We have observed a large sample of lo- 
cal Hipparcos clump giants at high spectral resolution and high 
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signal-to-noise ratio, and measured their metallicity and ele- 
mental abundances (Mishenina et al. 2006). Combined with a 
compilation of other studies providing metallicities of nearby 
clump giants, we have built a large unbiased sample of local 
giants to investigate the kinematical and chemical distributions 
of these stars. Our previous sample of distant giants was based 
on high resolution, low signal-to-noise spectra for 387 stars, 
spanning distances up to z=800 pc above the galactic plane, in 
the direction of the North Galactic Pole (NGP). The new distant 
sample now includes 523 stars up to z=l kpc, with improved 
distance and metallicity determinations. 

These two improved samples, local and distant, have also 
been used for other purposes, presented in separate papers. 
Kovtyukh et al. (2006) use the local sample to establish an 
accurate temperature scale for giants using line-depth ratios. 
Mishenina et al. (2006) investigate mixing processes in the 
atmosphere of clump giants. Finally Bienayme et al. (|2005), 
hereafter Paper III, apply two-parameter models on the com- 
bination of the local and distant samples to derive a realistic 
estimate of the total surface mass density within 0.8 kpc and 
1.1 kpc from the Galactic plane, respectively £o.8 kpc = 59-67 
M pc~ 2 and £ ukpc = 59-77 M©pc- 2 . 

Here we use these new data to focus on local properties 
of the thin disk which are so important to constrain its chemi- 
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cal and dynamical evolution : metallicity distribution, vertical 
metallicity gradient, age - metallicity relation (AMR) and age 
- velocity relation (AVR). Numerous studies of these proper- 
ties have been published, with however considerable disagree- 
ments reflecting the variety of tracers (open clusters, plane- 
tary nebulae, field dwarfs), discrepant metallicity scales, dif- 
ferent age determinations, or selection biases. A major contri- 
bution on the subject comes from the Geneva-Copenhagen sur- 
vey of the Solar neighbourhood by Nordstrom et al. (2004), 
which includes stellar parameters similar to ours, but for a 
much larger sample of dwarfs, and with photometric, less re- 
liable, metallicities. In the present work, the use of distant gi- 
ants allows us to probe larger distances above the galactic plane 
where kinematical distributions are no longer affected by lo- 
cal streams and moving groups, as studied by Famaey et al. 
(2005). Moreover, giants are well suited for age determina- 
tions, as shown in da Silva et al. (2006). We use their Bayesian 
method with isochrone fitting to compute ages and, similarly 
to them, we use the complete resulting probability distribution 
function of each star to bin the age axis. The combination of 
this pertinent method with the fact that we use spectrocopic 
metallicities for a large, homogeneous and complete sample, 
with well defined boundaries in magnitude and colour, should 
ensure that the new relations that we obtain are quite reliable. 
We have also computed for each star its probability of belong- 
ing, on kinematical criteria, to the thin disk, the thick disk, the 
Hercules stream and the halo, in order to reject the most prob- 
able non thin disk stars. 

Sections[2]and[3]describe the local and distant samples. We 
give details on the TGMET method and the new reference li- 
brary which have been used to improve the determination of 
J e ff, log g, [Fe/H], and M v for the distant giants observed at 



high spectral resolution, but low signal-to-noise (Section 3.2 1 



Sections [4] and [5] describe the determination of ages, Galactic 
orbits and population membership. Then we select the most 
probable thin disk clump giants and demonstrate the existence 



of a vertical metallicity gradient (Section 6.1 1. We present the 
AMR derived from the same stars in Section |6.2| while in 
Section |6.3| we discuss the AVR in U, V and W derived from 
a larger sample of clump giants where the most probable thick 
disk, Hercules stream and halo members have been rejected. 

2. The local sample of Hipparcos giants 

The sample of local giants, dominated by clump giants, con- 
sists of the 381 single Hipparcos field stars which follow the 
criteria : 

7r > 10 mas 
Sicks ^ -20° 
0.7 < B- V < 1.2 

My < 1.6 

where n is the Hipparcos parallax and Sjcrs the declination. 
It is thus a complete sample. The Johnson B-V colour has been 
obtained from the Tycho2 B T - Vj colour applying Eq. 1.3.20 
from ESA ( fT997l > : 

B - V - 0.850 (B T - V T ) 



Absolute magnitudes M v were computed with V appar- 
ent magnitudes resulting from the transformation of Hipparcos 
magnitudes H p to the Johnson system, calibrated by Harmanec 
(fT998l . 

Radial velocities have been mainly compiled from observa- 
tions on the ELODIE spectrograph at Observatoire de Haute- 
Provance (OHP). Some 177 local giants have been observed for 
this project (Mishenina et al. 2006), while radial velocities of 
other stars were retrieved from the ELODIE library (Prugniel 
& Soubiran |200T1 120041 and the ELODIE archive (Moultaka et 
al. 12004). For the remaining stars, we found radial velocities in 
Famaey et al. (2005 ) and Barbier-Bossat et al. (2000). In sum- 
mary, we have retrieved radial velocities for 220 stars in the 
various ELODIE datasets, for 54 stars in Famaey et al's cata- 
logue, for 107 stars in Barbier-Bossat et al's catalogue. We have 
also retrieved from these different sources information about 
the binarity of the stars. We have flagged 30 suspected spectro- 
scopic binaries presenting an enlarged or double peak of their 
cross-correlation function. 

Atmospheric parameters (r e ff, log g, [Fe/H]) have been 
compiled from the [Fe/H] catalogue (Cayrel de Strobel et al. 
2001 ) updated with a number of recent references. The [Fe/H] 
catalogue is a bibliographical compilation which lists deter- 
minations of atmospheric parameters relying on high resolu- 
tion, high signal-to-noise spectroscopic observations and pub- 
lished in the main astronomical journals. We have added to the 
compilation effective temperatures determined by Alonso et al. 
J200TT ). di Benedetto ( fT99"8l , Blackwell & Lynas-Gray JT9981 
and Ramirez & Melendez (2005). A number of other recent 
references providing spectroscopic (T^, log g, [Fe/H]) have 
been added to the [Fe/H] catalogue in an effort to keep it up 
to date. For the present work, the largest contributions come 
from Mishenina et al. (2006) for 177 stars and da Silva et al. 
(2006]) for 14 stars. For the older references, which were al- 
ready in Cayrel de Strobel et al. (2001), the largest contribu- 
tion comes from Mc William ( 1990) for 233 stars. This compi- 
lation provided [Fe/H] for 363 stars, adopting a weighted av- 
erage when several values where available for a given star (a 
higher weight was given to the most recent references). For 5 
remaining stars, an ELODIE spectrum was available, enabling 
the determination of atmospheric parameters with the TGMET 
method (see next section). We thus have just 13 stars which 
lack atmospheric parameters, representing 3% of the whole lo- 
cal sample. 

Combining atmospheric parameters from different sources 
can be a source of errors if some verifications are not made. 
Not all authors of spectroscopic analyses use the same temper- 
ature scales, Fe lines, and atomic data so that systematic differ- 
ences may occur in the resulting metallicities. In the present 
work, our narrow ranges in colour and luminosity suggest 
we deal with a very limited range of atmospheric parameters 
where temperature determinations from different methods usu- 
ally agree well. This is confirmed in our sample where 99 stars 
have at least two different determinations of r e ff. Computing 
the mean T e g for each of these 99 stars, we find standard de- 
viations ranging from to 140 K, with a median value of 40 
K, which is below the commonly admitted external error on ef- 
fective temperatures (~ 50-80 K). Only 6 stars show T eff deter- 
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minations deviating by more than 100 K. Similar verifications 
were made on [Fe/H] : the median value of standard deviations 
around the mean for stars having at least two determinations is 
0.09 dex. 

Hipparcos proper motions and parallaxes have been com- 
bined with radial velocities through the equations of Johnson 
& Soderblom (1987]) to compute the 3 velocity components 
(U, V, W) with respect to the Sun (the U axis points towards 
the Galactic Center). 

Figure [TJ shows the distribution of this sample in the planes 
M v vs T eB , M v vs [Fe/H] and V vs U. 

3. The distant NGP sample 

The distant sample has been drawn from the Tycho2 cata- 
logue (H0g et al. 2000). We have applied similar criteria as in 
Soubiran et al. (2003 ) to build the list of red clump candidates, 
just extending the limiting apparent magnitudes to fainter stars. 
A detailed description of the sample can be found in Paper III. 
The resulting sample consists of 523 different stars on a 720 
square degree field close to the NGP. The Tycho2 catalogue 
provides accurate proper motions and V magnitudes. High res- 
olution spectroscopic observations on ELODIE allowed us to 
measure radial velocities, spectroscopic distances and metal- 
licities. 

3.1. Spectroscopic observations, radial velocities 

The observations were carried out with the echelle spectro- 
graph ELODIE on the 1 .93 m-telescope at the Observatoire 
de Haute Provence. The performances of this instrument are 
described in Baranne et al. (1996). Compared to our previ- 
ous study in Paper I, 141 additional spectra have been ob- 
tained in Febuary and March 2003. The resulting 540 spec- 
tra cover the full range 390 - 680 nm at a resolving power 
of 42 000. The reduction has been made at the telescope with 
the on-line software which performs the spectrum extraction, 
wavelength calibration and measurement of radial velocities 
by cross-correlation with a numerical mask. The radial veloc- 
ity accuracy is better than lkms -1 for the considered stars 
(K stars). Our sample spans radial velocities from -139 to 85 
kms~' with a mean value of -12.8 kms~'. The mean S/N of 
the spectra at 550 nm is 22. Some 17 stars have been observed 
twice. For 13 stars, the correlation peak was enlarged or dou- 
ble, indicating the probable binarity of these stars which were 
flagged. 

3.2. Stellar parameters (T eS , log g, [Fe/H], M v ) 

We have performed the determination of stellar parameters T e s, 
log g, [Fe/H] and M v from ELODIE spectra using the code 
TGMET (Katz et al. [T9981 , like in Paper I. TGMET relies on 
the comparison by minimum distance of target spectra to a li- 
brary of stars with well known parameters, also observed with 
ELODIE (Soubiran et al. [19981 Prugniel & Soubiran |200T| ). 
As compared to Paper I, we have improved the content of the 
TGMET library because we were aware that the quality of 
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Fig. 1. Our local sample in the M v vs T e g, M v vs [Fe/H] and V 
vs U diagrams 
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TGMET results are very dependent of the quality of the em- 
pirical library which is used as reference. We present in this 
section the library that we built for the present study dealing 
with clump giants. We also present the tests which have been 
performed to assess the reliability of the TGMET parameters. 

The TGMET library must be built with reference spectra 
representative of the parameter space occupied by the target 
stars, with a coverage as dense as possible. The parameters of 
the reference spectra must be known as accurately as possible. 
Since our previous study of clump giants at the NGP, in pa- 
pers I and II, the TGMET library has been improved consider- 
ably. Many stars with well determined atmospheric parameters, 
compiled from the literature, and with accurate Hipparcos par- 
allaxes, have been added to the library as reference stars for 
T e «, log g, [Fe/H] and My. In particular the Hipparcos giants 
observed with ELODIE to build the local sample and analysed 
by Mishenina et al. (2006 ) have been added to the library. Fig. [2] 
shows the distribution of the TGMET library used for this study 
in the plane ([Fe/H], My). The clump area is densely covered 
down to [Fe/H] = -0.80. 

A small part of the TGMET library is presented in Table 
[T] The full Table is only available in electronic form, at the 
CDS. The calibrated Echelle spectra can be retreived from the 
ELODIE archived 



TGMET library 





1 

* 


I I 

■ * 






* * 
* ■ 
• 




■ 


. ■•*.•" 

• '.•r't 

• • " » 




-2 


-1 
[Fe/H] 



Fig. 2. Distribution of the TGMET library used in this study 
(724 reference stars observed with ELODIE) in the plane 
metallicity - absolute magnitude. 

In order to verify the TGMET results, we have run the code 
on ELODIE spectra of stars chosen in the library, with the best 
known parameters, degraded to a S/N typical of our target spec- 
tra (i.e. S/N=20). We have applied a bootstrap method : each 
test spectrum was removed in turn from the library, degraded 

1 http://atlas.obs-hp.fr/elodie/ 



to S/N=20, and its parameters determined by comparison to the 
rest of the library. To check results on My, we have selected the 
158 stars of the library with a relative error on their Hipparcos 
parallax lower than 10% and with 0.9 < B-V< 1.1. For [Fe/H] 
we have selected 199 stars with 0.9 < B-V< 1.1 having the 
most reliable spectroscopic metallicity determinations found in 
the literature. My and [Fe/H] determined from TGMET were 
then compared to their Hipparcos and literature counterpart, as 
shown in Figs. [3] and [4] The rms of the comparison, respec- 
tively 0.25 mag and 0.13 dex on My and [Fe/H], measure the 
accuracy of the TGMET results at S/N=20. The rms on My 
corresponds to an error in distance of 12%. 




(Hipparcos) 

Fig. 3. Comparison of TGMET absolute magnitudes from de- 
graded spectra to those deduced from Hipparcos parallaxes for 
a subset of 158 reference stars. 
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Fig. 4. Comparison of TGMET metallicities from degraded 
spectra to those from the literature for a subset of 199 refer- 
ence stars. 
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Table 1. Data and stellar parameters compiled for the TGMET library : HD/BD number, date of observation, T e g, log g, [Fe/H], 
My, quality flags qt, qf and qm for r e ff, [Fe/H] and My respectively (0 : data not available, 1 : poor, ... 4 : high), S/N of ELODIE 
spectrum at 550 nm, radial velocity, B — V, spectral type. The full table for all 724 stars is only available electronically at the 
CDS. The corresponding Echelle spectra can be retreived from the ELODIE archive. 
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In order to test the internal precision of TGMET on [Fe/H], 
we have compared the results obtained for the 17 stars ob- 
served twice (Fig. As can be seen, the agreement is very 
good (rms=0.05 dex). 



and deduced from the bootstrap test on degraded spectra, in 
0.25 mag bins. Similarly, Fig.|7]shows the two metallicity his- 
tograms, from the literature and from the bootstrap test. These 
histograms are perfectly aligned and present similar disper- 
sions which guarantees the lack of bias in the TGMET results. 



^0.8 -0.6 -0.4 -0.2 0.2 

[Fe/H] (lowest S/N) 

Fig. 5. Comparison of the TGMET metallicities obtained for 
the 17 target stars observed twice (rms=0.05 dex). 

An important verification has to be made to check that 
TGMET does not introduce a bias in the absolute magnitude 
and metallicity distributions of giants. In the following sec- 
tions, parameters of distant giants, relying on TGMET, will be 
compared to parameters of local giants, relying on literature 
and Hipparcos data. We thus have to ensure that these parame- 
ters are on the same scales. Fig.|6]shows the histograms of abso- 
lute magnitudes of the library's giants deduced from Hipparcos 




My 



Fig. 6. Absolute magnitude histograms of the library's giants 
deduced from Hipparcos (filled) and deduced from the boot- 
strap test on degraded spectra (red line). 
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Fig. 7. Metallicity histograms of the library's clump giants de- 
duced from the literature (filled) and deduced from the boot- 
strap test on degraded spectra (red line). 
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3.3. Distances, spatial velocities 

Distances have been computed for all the target stars from 
the TGMET My and Tycho2 Vj magnitude transformed into 
Johnson V. No correction of interstellar absorption was applied 
since it is supposed to be very low in the NGP direction. Proper 
motions, distances and radial velocities have been combined to 
compute the 3 velocity components (U, V, W) with respect to 
the Sun. 

Figure[8]shows the distribution of the 523 target stars in the 
planes M v vs r eff , M v vs [Fe/H] and V vs U. 



■ * t 



4. Ages, Galactic orbits 

Ages have been computed with the code PARAM developed by 
L. Girardi, available via an interactive web forrrj^] The method 
was initially developed by J0rgenson & Lindegren (2005) and 
slightly modified as described in da Silva et al. (20Q6J. It is a 
Bayesian estimation method which uses theoretical isochrones 
computed by Girardi et al. (2000) taking into account mass 
loss along the red giant branch. A convincing application of the 
method to derive the fundamental parameters of evolved stars 
in an open cluster is presented in Biazzo et al. (|2007). Inputs to 
be given to the code are the observed effective temperatures, ab- 
solute magnitudes, metallicities and related errors. The output 
for each star is a probability distribution function (PDF) of the 
age (and other parameters which are not used here). As shown 
in da Silva et al. (J2006), in their Fig. 5, the PDF of ages can 
be asymetric or even double peaked, especially in the case of 
red clump giants. As a consequence, ages are accurate for only 
a tiny part of our sample. This should be kept in mind for the 
use of individual ages. Nevertheless the ages have significance 
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Fig. 8. The NGP sample in the M v vs T eB , M v vs [Fe/H] and V 
vs U diagrams. 
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when used statistically. As a proof, the age-metallicity plot for 
the 891 stars (Fig. [9]l shows a regular trend and a remarkably 
low dispersion as compared to other studies (e.g. Nordstrom et 
al. 120041 da Silva et al. l2006l The 143 stars (83 local, 60 dis- 
tant) with relative age errors < 25% have been highlighted in 
Fig. [9] Considering only these stars, we measure a mean metal- 
licity of -0.06 with a dispersion of 0.10 dex for stars younger 
than 2 Gyr, whereas the mean metallicity of older stars (age > 
8 Gyr) is -0.44 with a dispersion of 0.27 dex. There is no young 
star with a metallicity lower than -0.32, and no old star with a 
metallicity higher than -0.13, contrary to common findings in 
samples of dwarfs, as for instance in Feltzing et al. (2001 ) and 
Nordstrom et al. d2004l > . It is important to note this property 
of our sample, because the existence of old metal-rich stars is 
often mentioned to explain the large dispersion of the AMR 
(Haywood 120061 1. We come back to the AMR of the thin disk 
in Sect.|6l 



Fig. 9. Age - metallicity diagram for the 891 stars. Stars (83 
local, 60 distant) with relative age errors lower than 25 % are 
highlighted as large filled circles. 
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The orbital parameters have been computed by integrat- 
ing the equations of motion in the galactic model of Allen & 
Santillan (1993), adopting a default value of 4 Gyr as the in- 
tegration time. The adopted velocity of the Sun with respect 
to the LSR is (9.7, 5.2, 6.7) kms 1 (Bienayme 1999), the so- 
lar galactocentric distance Rq = 8.5 kpc and circular velocity 
V LSR = 220 kms" 1 . 



5. Population membership 

The U vs V velocity distributions of the local and distant sam- 
ples can be compared from Figures [TJ and [8] It is clear, from 
these plots, that the two samples contain different kinematical 
populations. In the local sample, the velocities are clumpy 
and reflect moving groups and superclusters that dominate the 
kinematics in the solar neighbourhood. Compared to Fig. 9 of 
Famaey et al. (2005), we can identify the Hercules stream at 
(U, V) - (-40, -50) kms -1 , the Hyades-Pleiades supercluster 
at (U, V) =s (-30,-20) kms -1 , and the Sirius moving group 



at ( U, V) (0, 0)kms . There are very few high velocity 
stars that could correspond to the thick disk. On the contrary, 
the velocities of the distant sample are better mixed with 
higher dispersions. This reflects the dynamical heating of the 
disk together with the growing number of thick stars with 
increasing distance to the plane. In order to build a sample of 
pure thin disk stars, we have performed the classification of all 
the stars into different kinematical populations. We have taken 
into account the Hercules stream because its velocity ellipsoid 
is just intermediate between that of the thin disk and the thick 
disk, and is likely to contaminate both populations. We did not 
attempt to distinguish the other groups of the thin disk. 

We assign to each star its probability to belong to the thin 
disk, the thick disk, the Hercules stream and the halo on the ba- 
sis of its (U, V, W) velocity and the velocity ellipsoids of these 
populations, in the same way as Soubiran & Girard (2005) and 
with similar kinematical parameters of the populations. In the 
distant sample we find that 305 stars and 65 stars have a prob- 
ability higher than 80% to belong to the thin disk and the thick 
disk respectively. In the local sample, the numbers are 304 and 
11. 

One important question that we can immediately investi- 
gate thanks to this kinematical classification is whether the thin 
disk and the thick disk overlap in age and metallicity. Our data 



strongly suggest that this is the case. Fig. 10 shows with differ 



ent symbols the age-metallicity diagram for the most probable 
thin disk and the thick disk stars, restricted to relative age er- 
rors lower than 25 % (suspected binaries rejected). It is clear 
that the oldest thin disk stars and thick disk stars overlap in the 
metallicity range -0.30 < [Fe/H] < -0.70, and age range 8-10 
Gyr. It is also worth noticing that there are no young thick disk 
stars. 



Fig. 10. Age - metallicity diagram for stars with relative age 
errors lower than 25 % and belonging to the thin disk (crosses) 
and the thick disk (filled circles). 




All the parameters which have been determined as de- 
scribed in the previous sections are presented in Table [2] The 
full table with all 891 stars is only available in electronic form 
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at the CDS. The file with the age PDFs is also available upon 
request. 
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6. The thin disk traced by clump giants 

Among the many studies that can be done with the new sam- 
ple presented here, we focus on the properties of the thin disk 
probed for the first time up to large distances above the Galactic 
plane, from a complete stellar sample and with 3D kinemat- 
ics and spectroscopic metallicities. In order to work with a ho- 
mogeneous sample, with well defined boundaries in both its 
local and distant counterparts, we have selected clump giants 
on the basis of a colour and absolute magnitude restriction : 
0.9 < B -V < 1.1, 0.0 < M v < 1.6. According to Koen 
& Lombard (2003), this ensures the lowest contamination by 
other giants. Rejecting suspected binaries, 597 stars fall into 
these limits. We further restrict the sample to the 396 stars hav- 
ing a probability higher than 80% to belong to the thin disk. 
In this section we investigate some basic distributions of this 
sample. 



6.1. Raw metallicity distribution and vertical gradient 

We compare the metallicity distributions of the local and dis- 
tant clump giants in Fig. 1 1 The local sample has an average 
of [Fe/H]=-0.11 and a standard deviation of <T[F e /H] = 0.15 



whereas the distant sample has an average of [Fe/H]=-0.21 and 
a standard deviation of cr^/^ = 0.17. The metallicity dis- 
tribution of the thin disk is thus significantly shifted towards 
lower values at larger distance above the galactic plane. This 
is not due to the comparison of metallicities coming from the 
literature for the local sample and from TGMET for the distant 
sample since we have verified that the two scales are consis- 



tent (Sect. 3.2 1. More likely this difference indicates a vertical 



metallicity gradient which is represented in Fig. 12 using as the 



distance the maximum height from the plane, Zmax, reached 
by the star in its galactic orbit. A linear fit indicates a gradient 
of 3[Fe/H]/3Z = -0.31 ± 0.03 dex kpc 1 . Taking for each star 
its current distance from the plane, instead of Zmax, leads to a 
consistent result of <9[Fe/H]3z = -0.30 + 0.03 dex kpc 1 . 

According to numerous previous studies, the existence of a 
vertical metallicity gradient in the thin disk seems to be firmly 
established. However the value of its amplitude, constrained by 
the observation of different kinds of tracers, at various distances 
from the Sun, is still oscillating between ~ -0.25 and -0.35 dex 
kpc~'. Using open clusters, Piatti et al. (Q995) find -0.34 dex 
kpc -1 whereas Carraro et al. (1998) measure -0.25 dex kpc 1 
and Chen et al. J20031 measure -0.295+0.050 dex kpc 1 . Other 
studies are based, like ours, on field stars and have used kine- 
matical information to select thin disc stars. This is the case of 
Marsakov & Borkova (2006 ) who have selected the most prob- 
able thin disk stars in their compilation of spectroscopic abun- 
dances, using their 3D velocties and orbital parameters. They 
measure a gradient of -0.29 + 0.06 dex kpc -1 . Bartasiute et al. 
(2003) have observed 650 stars at high galactic latitude, up to 
1 . 1 kpc, and identified thin and thick disk stars on the basis of 
their rotational lag. They measure a gradient of -0.23 + 0.04 
dex kpc 1 in the thin disk. 

The direct comparison of the metallicity distribution of 
our sample with other distributions probing different galac- 
tic volumes would imply a scaleheight correction. The reason 



is that metal-poor stars, which have hotter kinematics, have 
a larger scaleheight than more metal-rich stars, and may be 
under-represented in local samples. A correction, relying on 
a mass model of the disk, would thus increase the number of 
metal-poor stars with hotter kinematics which would have been 
missed in our sample. On the contrary, stars more metal-rich 
than the Sun are supposed to be over-represented in local sam- 
ples (see for instance Fig. 3 in Haywood 2006). We have not 
attempted to correct such bias in our sample and we restrict 
here the discussion to a qualitative comparison between dwarfs 
and giants. 



When we compare the metallicity distribution of clump gi- 
ants to that of dwarfs, as presented by Haywood (2002}, we 
find a good agreement for the metal-poor side. We confirm 
with this new sample Haywood's finding that the thin disk is 
not an important contributor to stars with [Fe/H] < -0.5. We 
find that 2.5% of our sample has [Fe/H] < -0.5 with the most 
metal-poor thin disk giant at [Fe/H] =-0.71. According to Fig. 
3 in Haywood (2006 ), the scaleheight correction factor is com- 
prised between ~ 1.5 and 3.5 in the metallicity range -0.70 < 
[Fe/H] < -0.50. Taking this correction into account would not 
change substantially our findings. 



On the contrary, we find a significant difference between 
clump giants and dwarfs for the metal-rich side of the [Fe/H] 
histogram. Haywood (2002) finds that 40-50% of long-lived 
dwarfs have a metallicity higher than [Fe/H]=0 whereas the 
proportion is only 20% in our local sample and 13% in 
our distant sample. Super Metal-Rich ([Fe/H] > +0.20) FGK 
dwarfs are quite usual in the Solar Neighbourhood whereas 
we have only two thin disk clump giants at [Fe/H] =+0.21 and 
[Fe/H] =+0.27. Our first guess was that such a low ratio of 
metal-poor stars and metal-rich stars was correlated with the 
colour cuts that we used to restrict the sample to clump gi- 
ants. We have verified that this is not the case by comparing 
the metallicity histograms of local giants (0 < M v < 1 .6) in 
the B-V intervals [0.9; 1.1] and [0.7; 1.2]. We found that the 
metal-poor sides are similar. The ratio of metal-rich stars turns 
out to be slightly higher in the extended colour interval : 24% 
instead of 20%. We conclude that our adopted colour cutoff af- 
fects the metallicity distribution in a way that metal-rich stars 
are slightly under-represented. This bias is however not suffi- 
cient to reconcile the metallicity distribution of clump giants 
with that of dwarfs. 



Pasquini et al. (2007) have also noticed a difference in the 
metallicity distribution of giants and dwarfs hosting planets. 
They propose as an interpretation the pollution of stellar atmo- 
spheres, causing a metal excess visible in the thin atmosphere 
of dwarfs, while diluted in the extended envelope of giants. Our 
sample suggests that the difference is not limited to stars host- 
ing planets so that the pollution hypothesis should be inves- 
tigated in a more general context. If validated in the general 
case, it would imply that dwarfs are not appropriate to probe 
the chemical history of the Galaxy. 
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Fig. 11. Metallicity distribution of thin disk clump giants of the 
local (filled) and distant samples (red line). 




0.5 



[Fe/H] 



Fig. 12. Vertical gradient in the metallicity distribution of thin 
disk clump giants. 




6.2. Age - metallicity relation 



We use the same method as da Silva et al. (2006) to determine 
the AMR of the sample. For each time interval (A,=l Gyr), we 
measure the cumulative [Fe/H] by adding the measured [Fe/H] 
of each star weighted by its PDF. With such a method, a given 
star contributes to several bins which are consequently not in- 
dependant. However it is a good way to take errors on ages into 
account. The mean metallicity and dispersion per age bin are 
presented in Fig. 13 and Table[3] A remarkable result is the low 
dispersion obtained at all ages. Subtracting the estimated ob- 
servational error (0.09 dex for local stars, 0.13 dex for distant 
stars) yields to a cosmic scatter in [Fe/H] lower than 0.12 dex. 
A transition occurs around 4 Gyr in both the mean metallicity 
and dispersion. From 10 Gyr to 4 Gyr, we see a very smooth 
and regular increase of the mean metallicity, 0.01 dex per Gyr, 
with constant spread, which characterizes an homogeneous in- 



terstellar medium. An upturn occurs at 4 Gyr with a steeper 
metallicity rise at younger ages. 

What is the state-of-the-art of the AMR in the Solar 
Neighbourhood and how does our relation compare with previ- 
ous ones ? Despite numerous studies on the subject over nearly 
30 years, there is still no consensus on the existence or not 
of a slope in the AMR, neither in the amplitude of the cos- 
mic scatter. Major contributions like Edvardsson et al. (JT993 ), 
Feltzing et al. (2001) and Nordstrom et al. (2004), using classi- 
cal isochrone ages, find little evidence for a slope in the relation 
of [Fe/H] with age, and a broad dispersion (<X[p e / H ] > 0.20 dex). 
In contrast, Rocha-Pinto et al. (2000), using chromospheric 
ages, find a significant trend in the AMR, with lower disper- 
sion (o~ [Fe/H] ~ 0.12 dex). Pont & Eyer (2004) have re-analysed 
Edvardsson et al.'s sample with a Bayesian approach and also 
find a significant trend with a dispersion cr [Fej / H j < 0.15 dex. 
We note that all these studies involve nearby dwarf stars. To 
our knowledge, the only AMR based on giants is that of da 
Silva et al. (2006). Despite the small size of their sample, they 
find like us metallicities rising from [Fe/H] ~ -0.23 at 10.5 Gyr 
to [Fe/H] ~ 0.00 at 0.5 Gyr. The shape of their AMR is how- 
ever different of ours, shallower at young ages and steeper at 
old ages. The dispersion of their AMR is also much larger than 
ours, reaching 0.30 dex in the oldest age bins. 

We notice that the rather large metallicity variation that we 
observe in the 4 youngest bins in our AMR is also visible in 
the AMR derived by Nordstrom et al. (2004) and by Feltzing et 
al. ([2001). Both studies interpret this feature as a bias against 
young metal-poor dwarfs due to a colour cut. This explana- 
tion is not valid for our sample since we have verified that our 
colour cuts only affect very slightly the metal-rich part of the 
metallicity distribution (see previous Sect.). We thus conclude 
that this peculiar shape of the AMR is real. 

Piatti et al. (119951) and Carraro et al. 0998) have corrected 
their AMR from the positional dependency, justified by the use 
of open clusters. Open clusters have a wide spatial distribution 
and trace different histories of the chemical evolution, depend- 
ing on their galactocentric distances. The AMR has thus to be 
corrected from the observed radial metallicity gradient, which 
has an amplitude of 0.07 dex kpc~' according to Piatti et al. 
dl9951 l. or 0.09 dex kpc 1 according to Carraro et al. (11998b . 
Field stars are also supposed to be affected by a radial metallic- 
ity gradient. A consequence of orbital diffusion is that samples 
of nearby stars may include stars born in the inner or outer parts 
of the disk where the chemical enrichement may have been dif- 
ferent from that of the Solar Neighbourhood. Such stars are eas- 
ily recognized with their orbital parameters Rmin and Rmax, 
respectively perigalactic radius and apogalactic radius, differ- 
ent from that of the true local stars. It is worth mentioning that 
Edvardsson et al. (1993) have studied the AMR for stars re- 
stricted to the solar circle and still found a large and significant 
scatter. Our sample of thin disk clump giants is free from the 
influence of stars from other galactocentric distances since our 
kinematical selection has naturally eliminated stars on eccen- 
tric orbits. 

The question whether the AMR should be corrected from 
the vertical metallicity gradient is more difficult to assess. We 
note that Carraro et al. (1998]) have not attempted to correct 
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Table 3. Age-metallicity relation derived from our sample of 
thin disk clump giants. N is the number of stars contributing 
to each age bin. It is fractional because we use the complete 
probability function of each star to bin the age axis (see text). 



spanning larger distances from the Galactic plane, is well suited 
to investigate these relations. However, for such a purpose, we 
cannot work on the thin disk sample which was built to study 
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their open cluster AMR from the observed vertical metallicity 
gradient. Moreover they argue that "In the case of field stars, 
orbital diffusion is expected to be effective enough to smooth 
out a vertical metallicity gradient within a single-age popula- 
tion, so that the vertical structure of the disk is dominated by 
the different scaleheights of different age populations". We also 
note that, in the case of field star AMRs, while the radial mi- 
gration is often refered to (Edvardsson et al. 119931 Haywood 
2006 ), the influence of the vertical metallicity gradient is not 
discussed. 



Fig. 13. Age - metallicity relation of thin disk clump giants. 
The error bars represent the dispersion in each bin, including 
observational errors and cosmic scatter. 




age (Gyr) 



6.3. Age - velocity relation 

The thin disk AVR has been revisited recently by Seabroke & 
Gilmore (120071) using the data of Nordstrom et al. (|2004 ) and 
Famaey et al. (2005 ). They show that the kinematical streams 
in these local samples do not permit one to safely constrain the 
relations in the U and V directions, contrary to the W direction 
where the samples are well mixed. Our sample of clump giants, 



the metallicity and age distributions, in Sections 6.1 and 6.2 
Our selection of thin disk stars on a kinematical criterion has 
favoured stars in the central parts of the velocity ellipsoid, with 
moderate velocities, resulting in a serious kinematical bias. In 
order to study how the velocity dispersions increase with time, 
we need to work also with the warmer part of the thin disk, 
but excluding as well as possible stars which do not follow the 
kinematical behaviour of the thin disk. To do so, we consider 
our distant sample of clump giants and reject stars having a 
probability higher than 80% to belong to the thick disk, the 
Hercules stream and the halo, resulting in 320 stars. Results 
are presented in Fig. 14 and Table[4] 

An important question is whether the dispersions saturate 
at a given age of the thin disk. Seabroke & Gilmore (2007) 
have shown that local data are in agreement with several mod- 
els of disk heating : continuous or with saturation at 4.5, 5.5 
and 6.5 Gyr. Our data show evidence for a transition at ~ 5 
Gyr, with saturation occuring in V at 29 kms -1 and in W at 
24 kms -1 . The velocity dispersion in U seems to increase 
smoothly, reaching 46 km s at 9.5 Gyr. A consequence is that 
the velocity ellipsoid axis ratios o-y/cru and (Twjcru are not 
constant. The ratio <TvI(Tu is related to the Oort constants and 
is expected to be ~0.5. Here this ratio is varying from 0.55 at 1- 
2 Gyr to a maximum value of 0.68 at 4-5 Gyr. The ratio cr w /cru 
is related to the scattering process responsible of the dynami- 
cal heating of the disk. With our data, it has a maximum value 
of 0.56 at 4-5 Gyr. Although these ratios are supposed to be 
constant in an axisymmetric Galaxy, there are previous reports 
of variations related to colour or spectral types (e.g. Mignard 
120001) . 

We mention the study by Vallenari et al. (|2006 ) who have 
also probed the thin disk kinematics towards the NGR Their 
method is however significantly different from ours since they 
analyse, through a galactic model, proper motions and the 
colour magnitude diagram of ~ 15000 stars down to V=20. 
Their best-fit for the velocity dispersions of the thin disk, pre- 
sented in 4 age bins, differ significantly from ours, especially 
in the oldest age bin (7-10 Gyr) where their values are lower by 
3cr. 

Simple statistics on our sample gives (cry, o>, cr w ) = 
(41.5, 26.4, 22. l)kms _I , significantly higher than values de- 



termined from late-type Hipparcos stars (e.g. Bienayme 1999, 
Mignard 120001 ). Although we cannot rule out the contamina- 
tion of the sample with thick disk stars, it nicely compares to 
recent results by de Souza & Teixeira (2007 ) who show that 
Mignard's sample is better explained by the superposition of 
2 velocity ellipsoids the hotter one having (ov, cry, crw) = 
(41.0,27.0, 19.0) km s -1 . It is also worth noticing in Table 
that the mean U and W are roughly constant at all ages whereas 
V declines from — 14 kms~' to -21 kms -1 . We retrieve for 
U and V the Solar motion with respect to late-type stars, as 
determined by Mignard (2000), although we find a significant 
difference in W. We get a mean value of W = 11.5 kms~' , 
whereas he finds values around 7 km s . We recall that our W 
velocities of the distant stars at the NGP rely mainly on radial 
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Table 4. Age-velocity relation derived from 320 distant clump 
giants, the most probable thick disk, Hercules stream and halo 
stars being excluded. 
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velocities, which have an accuracy better than 1 km s _1 , and 
thus are not affected by uncertainties on distances and proper 
motions. 



Fig. 14. Age - velocity relations of distant clump giants, the 
most probable thick disk, Hercules stream and halo stars being 
excluded. 
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7. Summary 

The data presented here are the result of several years of ef- 
fort to obtain high resolution spectra for a large and complete 
sample of clump giants. Besides our own observations, on the 
ELODIE spectrograph at OHP, we have also taken advantage 
of other available material like the Hipparcos and Tycho2 cata- 
logues, the [Fe/H] catalogue (Cayrel de Strobel et al. 2001) up- 
dated with a number of new references and the PARAM code to 
derive ages (da Silva et al. l20061 >. We have described how these 
data were combined to provide a catalogue of stellar parame- 
ters for 891 stars, mainly giants, giving atmospheric parameters 
with spectroscopic metallicities, absolute magnitudes and dis- 



tances, galactic velocities (U, V, W), orbits, ages and population 
membership probabilities. 

Our main motivation in conducting this project was to 
probe the Galactic disk using an unbiased and significant sam- 
ple, with high quality data, in particular with spectroscopic 
metallicities and accurate distances and radial velocities. We 
have chosen to observe giants in the direction of the NGP in or- 
der to reach distances to the galactic plane, up to 1 kpc, which 
are not covered by spectroscopic surveys, usually limited to 
the closer Solar Neighbourhood. Clump giants are particularly 
well suited for this purpose. Compared to previous studies on 
the subject, our analysis presents several improvements, which 
are briefly outlined: 

- for binning the age axis, we have considered for each star its 
entire age PDF, instead of averaging it, following da Silva 
et al. (120061 

- we have considered several kinematical populations likely 
to be present in our sample: the thin disk, the thick disk, the 
Hercules stream and the halo 

- in order to study the thin disk metallicity and age distribu- 
tions, we have taken care to select stars with the highest 
probability of belonging to this population 

- in order to study the thin disk velocity distribution, we have 
taken care to reject the most probable non thin disk stars 

Our results are summarized as follows : 

- we do not find any young metal-poor stars nor old metal- 
rich stars, contrary to common findings in dwarf samples 

- the old thin disk and the thick disk overlap in the metallicity 
range -0.70 < [Fe/H] < -0.30 and age range 8-10 Gyr 

- among stars with accurate individual ages, we do not find 
any young thick disk stars 

- the metallicity distribution of our sample of thin disk clump 
giants extends down to [Fe/H]^-0.70, but the fraction of 
stars with [Fe/H]<-0.50 is only 2.5% 

- the metallicity distributions of giants and dwarfs differ sig- 
nificantly on the metal-rich side: metal-rich giants are less 
frequent 

- a vertical metallicity gradient is measured in the thin disk: 
<9[Fe/H]/3Z = -0.31 + 0.03 dex kpc 1 

- the AMR of the thin disk presents a low dispersion, imply- 
ing a cosmic scatter lower than 0.12 dex, in agreement with 
previous findings by Rocha-Pinto et al. (2000) and Pont & 
Ever d2"004l l 

- 2 regimes are visible in the AMR of the thin disk : from 
10 Gyr to 4 Gyr, the metallicity increases smoothly by 0.01 
dex per Gyr, while for younger stars the rise of [Fe/H] is 
steeper 

- in the thin disk, the V and W dispersions saturate at 29 and 
24 km s respectively at ~ 4-5 Gyr, whereas U shows con- 
tinuous heating 

- the Solar motion is found to be nearly constant in U and 
W with respect to stars of all ages, while the amplitude of 
the asymmetric drift increases from 14 to 21 kms -1 with 
respect to young and old stars respectively 

Acknowledgements. We are grateful to L. Girardi for computing 
the ages for the 891 stars of this sample. This research has made 
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